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1. Introduction 

Continuous time random walks (CTRWs) are random walks with random waiting 
times Wk between jumps </&. They have been applied in physics to a variety of systems 
exhibiting “anomalous diffusion,” with heavy-tailed waiting times leading to subdiffu- 
sive processes whose variance grows oc C 3 , 0 < f3 < 1, and with heavy-tailed jumps 
leading to superdiffusive processes which exhibit a faster scaling than Brownian motion 
(Metzler and Klafter, 2000). For a variety of applications, see e.g. Berkowitz et al. 
(2006); Henry and Wearne (2000); Fedotov and Iomin (2007); Raberto et al. (2002); 
Schumer et al. (2003). Scaling limits of CTRWs are non-Markovian time-changes 
of M d -valued Markov processes (Meerschaert and Scheffler, 2004; Kolokoltsov, 2009; 
Kobayashi, 2010). 

The main tool for the analysis and computation of the distribution of CTRW limits 
is the (fractional) Fokker-Planck equation (FPE; considered here as a synonym with 
Kolmogorov Forward Equation) (Barkai et ah, 2000; Langlands and Henry, 2005); for 
textbooks with a quick introduction to fractional derivatives see e.g. Meerschaert and 
Sikorskii (2011) or Kolokoltsov (2011). Governing FPEs have been derived in the 
literature, whose solutions can be roughly classified as follows: 

• Classical (strong) solutions (Kolokoltsov, 2009; Hahn et ah, 2010; Kolokoltsov, 
2011; Magdziarz et ah, 2014; Nane and Ni, 2015), where the derivation assumes 
that the underlying space-time Feller process (see below) has (differentiable) 
probability densities; 

• Solutions in Banach space settings, in the framework of fractional evolution 
equations (Priiss, 2012; Bajlekova, 2001; Baeumer and Meerschaert, 2001; Baeumer 
et ah, 2005; Umarov, 2015) where the derivation is based on the assumption that 
the coefficients do not vary in time; 
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• Mild solutions based on Fourier-Laplace transforms (Becker-Kern et al., 2004; 
Meerschaert and Scheffler, 2008; Jurlewicz et al., 2012). These allow for a cou¬ 
pling between jumps and waiting times but assume constant coefficients. 

One aim of this paper is to unify the above results and to derive governing FPEs without 
these restricting assumptions. 

A further important analytical and computational tool for anomalous diffusion pro¬ 
cesses is the (fractional) Kolmogorov backward equation. It may be used to calculate 
distributions of occupation times and first passage times for anomalous diffusion pro¬ 
cesses (Carmi et ah, 2010). In groundwater hydrology, scaling limits of CTRWs model 
the spread of contaminants in an aquifer (Berkowitz et ah, 2006; Schumer et ah, 2003), 
and (non-fractional) Kolmogorov backward equations have already been used to model 
the distribution of pollutant sources and travel times (Neupauer and Wilson, 1999). A 
mathematical framework for CTRW scaling limits and fractional Kolmogorov backward 
equations would hence be applicable to problems in groundwater hydrology, but has 
yet to be established, which is the second aim of this paper. 

The following topics are not discussed in this article in order to maintain the focus 
on governing equations, but they should be mentioned as they are closely related: 

• If each waiting time and jump pair (H4, Jk ) is coupled, their order is important: 

CTRWs assume that W *. precedes whereas OCTRWs (overshooting CTRWs) 

assume that J*, precedes Wk■ The scaling limits of these two processes may be 
as different as having mutually disjoint supports for all t > 0 (Jurlewicz et ah, 
2012; Straka and Henry, 2011). This paper focuses on limits of CTRWs. 

• In our analysis, we do not assume that any stochastic process admits a Lebesgue 
density, hence the FPE (Th 5.2) is given on the Banach space of positive mea¬ 
sures; such a result is apparently new. If the Feller process (A r , D r ) below admits 
a (suitably regular) density, then the CTRW limit does so, too (Magdziarz et ah, 
2014). In general, however, and in particular for the three examples discussed 
in the last section, the existence of densities is unconfirmed. 

• Our analysis defines CTRW limits via a continuous mapping approach (Theo¬ 
rem 2.1), and the underlying assumption is the convergence of Feller jump pro¬ 
cesses to a Feller diffusion process with jumps. If the sequence of Feller jump 
processes is specified, then the sequence of CTRWs is also specified, which can 
be illuminating for applications and the simulation of sample paths. We skip 
this content with a warning that convergence can be difficult to establish ( Jacod 
and Shiryaev, 2002; Kolokoltsov, 2011). 

Organization of this paper: In Section 2 below, CTRW scaling limits are introduced 
in a very general setting. Section 3 introduces the Banach space setting needed for 
the derivation of Kolmogorov backward equation (Section 4) and the FokkerAPlanck 
Equation (FPE, Section 5). Finally, Section 6 contains three examples from statistical 
physics which illustrate the forward and backward governing equations. 

2. Scaling limits of CTRWs 

We introduce CTRW limit processes by closely following Meerschaert and Straka 
(2014): Let c > 0 be a scaling parameter, and write A c (n) for the position after the 
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n-th jump, and D c (n) for the time of the n-th jump. We assume that after each jump, a 
CTRW is renewed. More precisely, ( A c {n + 1), D c {n + 1)) depends on the previous tra¬ 
jectory (kL c (0), ZU(0)),..., (A c (n), D c [n)) only through the latest pair (H c (n), _D c (n)); 
but this is equivalent to {(A c (n), D c (n))} n£ ^ 0 being a Markov chain with state space 
M'm. yy e assume that the sequence D c {n) is strictly increasing. 

By setting A c (t ) = H c (|_fJ), D c (t) = D c ([t\), a Markov chain as above defines a 
trajectory [0, oo) 3 i 4 (A c (t), D c (t)) £ M d+1 . This trajectory can then be mapped 
to a CTRW trajectory as follows: Define the right-continuous inverse E c (t) := inf{n : 
D c (u ) > t} of D c . Write A c _ for the left-continuous version of A c . Then the CTRW 
trajectory is given by 

X c (t) =A c _oE c _(t+), 

that is, by the right-continuous version of the composition of the two left-continuous 
processes A c _ o E°_ (Straka and Henry, 2011, Lemma 3.5). One may then exploit the 
Skorokhod continuity of this path mapping to obtain the CTRW scaling limit as c —> oo: 

Theorem 2.1. Suppose that as c —>■ oo, we have the weak convergence 

(2-1) {(VacrJ),B o (LcrJ))}^ 0 ^{(X.A.)} r > 0 

in the J\ topology on cadlag paths in R d+1 , where D r is a.s. strictly increasing and 
unbounded. Then we also have the weak convergence 

{X c (t)} tgR =t {W(f)} tgR 

in the J\ topology on cadlag paths in where 

(2.2) X(t) = A_ oE(t+), 

A _ denotes the left-continuous process {^4(t—) }t>o and E(t) = inf {;u : D u > t}. 

Proof. This theorem is a direct consequence of Proposition 2.3 in Straka and Henry 

( 2011 ). □ 

We stress that Elf) and X(t) are in general not Markovian. 

Due to the above theorem, the large class of possible CTRW limit processes is hence 
essentially given by (2.2) and an M. d+1 valued process (A r , D r ) which is the weak limit 
of a sequence of (continuous time) Markov chains, where D r is strictly increasing and 
unbounded. Such processes contain the class of diffusion processes with jumps, in the 
sense of Jacod and Shiryaev (2002). Details on the convergence of Feller-jump processes 
to a Feller diffusion process with jumps as in (2.1) are e.g. in Theorem IX.4.8 of the 
mentioned textbook, and in Kolokoltsov (2011) with somewhat more specificity. 

The idea that CTRWs are essentially random walks in space-time was seemingly first 
introduced explicitly to CTRWs by Weron and Magdziarz (2008), and used in Henry 
et al. (2010) to derive a Fractional Fokker-Planck Equation (FPE) with space- and 
time-dependent drift. (For a more detailed derivation of the FPE, see Magdziarz et al. 
(2014).) 

The following scaling limits ( A r , D r ) have been considered in the literature: Uncou¬ 
pled and coupled stable limits (Meerschaert and SchefHer, 2004; Becker-Kern et ah, 
2004), triangular array limits (Meerschaert and SchefHer, 2008; Jurlewicz et ah, 2012), 
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position-dependent, stable-like limits (Kolokoltsov, 2009) and stochastic differential 
equations with diffusion component A r and subordinator D r (Weron and Magdziarz, 
2008; Magdziarz et ah, 2014). 

To specify the class of space-time limit processes (A r , D r ), we first define the operator 
A 0 : Co(M d+1 ) -G C 0 (M d+1 ) (with Einstein notation) by 

(2.3) Aof(x, s ) = 6*(x, s)d Xi f(x, s ) + y(x, s)dj(x , s) + ^a tJ (x , s)d Xi d Xj f(x, s ) 

+ J j [f(x + z,s + w)^f(x,s)-z % l(\\z\\<l)d Xi f(x,t)]K(x,s-,dz,dw). 

2 gRd U!>0 

We adopt the following basic conditions on the coefficients are: for i,j = 1 ,...,d 
the mappings (x,s) H > b l (x,s), (x,s) H y a^(x,s), (x,s) H > y(x, s), are in Cfe(M d+1 ), 
the measures K(x,s]-,-) are Levy measures for every (x, s) G M d+1 and Kg(x,s) : = 
ff K(x, s; dz , dw)g(z, w) lies in Cfe(M d+1 ) for every g G i?b(M d+1 ) (bounded measurable) 
which is 0 in a neighbourhood of the origin (Jacod and Shiryaev, 2002). We note 
however that these conditions are not sufficient for Aq to generate a Feller process; for 
sufficient conditions, consult e.g. Applebaum (2009, Ch 6). 

We assume that (A r , D r ) is a Feller process with strongly continuous semigroup 
(T r ,r > 0) acting on C'o(M cZ+1 ). The infinitesimal generator A of (T r ,r > 0) is such 
that C'o(M d+1 ) C Dom(M) and Af = Aof for all / G C'o(M d+1 ); for details, see e.g. 
Ch 6.7 in Applebaum (2009). We write P x,s for the (canonical) probability measure 
induced by ( T r ,r > 0) and P x,s (Ao = x, Dq = s) = 1. The requirement that D r 
be strictly increasing a.s. means that y(x, s) > 0, that the diffusive component of 
D r vanishes, that the measures K(x,s]-,-) are supported on x [0, oo), and that 
Jq 1 wK(x, s, M d , dw) < oo. Moreover, the truncation term in the integral does not apply 
to the d + 1st coordinate. For technical reasons, we require another, not very restrictive 
assumption: 

Transience in the time-component: If / G C^R^" 1 " 1 ) has support Supp(/) C 
R d x (—oo, B] for some 5 gR, then the potential of /, 

/•OO 

(2.4) Uf(x,s)= T r f (x, s) dr 

Jo 

is a continuous function with U /(•, s ) G Cc^R 0 ') for all s; i.e., with a slight abuse 
of notation there exists a kernel U such that 

Uf(x,s)= j U(x,s]dz,dw)f(z,w). 

For example, if D r is a subordinator then this assumption is satisfied (Bertoin, 1999). 
U is commonly referred to as the potential kernel of the semigroup (T r ,r > 0). 

We can now give a result which characterises the distribution of X t for Lebesgue- 
almost every t G R: 
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Theorem 2.2. Let H(x,s;v) := K(x, s; M d , (y , oo)), v > 0, and assume the following 
uniform integrahility condition: 


sup H(x,s;v) ) dv < oo. 
lo V (s,s)eR d+1 J 

Moreover, for h(x,s ) = f(x)g(s) with f G Co(K d ) and g G C c (M) (compact support) 
define the linear maps 

(2.5) s) := h(x, s)-f(x, s) + / h(x, s + v)H(x, s;v) dv. 


' v>0 


(2.6) Th(x, s) := h(x, s)'j(x, s) + 


u>0 -/zgR d 

T/ien i/ie CTRW limit process X t from (2.2) satisfies 


h{x + 2 , s + v)K(x, s; cfc x (u, oo)) dv 


(2.7) [ E*’ s [f(X t )g(t)}dt = UVh(x,s). 

J t>S 

and the OCTRW limit 


E*’ s [f(Y t )g(t)}dt = Urh(x,s). 


f t>S 


Proof. First note that H(x,s]v ) is decreasing to zero on v G (0, oo), for every (x,s) G 
R d+1 , since it is the tail function of a Levy measure. Hence ^h(x, s) —» 0 as x, s —» ±oo. 
Furthermore d'h is continuous by the Dominated Convergence Theorem and its support 
bounded above in s. Hence U^/h is well defined. 

Let h(x,t) = f(x)g(t) for some non-negative / G Co(M d ) and g G C c (R). Then by 
Tonnelli’s theorem, continuity of Lebesgue measure and the jumps of X t being countable 
the left-hand side of (2.7) equals 


hgR 


g(t)E x ’ s [f(X t )]dt = E x ’ s 


9(t)f(X t ) dt 


_Jt£ R 




g(t)f{X t _) dt 


_Jt£ R 




11£ R 


Now multiply the equation in Theorem 2.3 of Meerschaert and Straka (2014) by g(t) 
(neglecting Y t , V t and R t ) and integrate over t G 1, to get 


<7(0E 1 -* [/(A',-)] dt 


he r 



f(y)l(y,t)U(x,s;dy,dt ) 


+ 



U(x, s; dy, dr)H(y, r; t - r)f(y) dt 


Note that we may replace u x,T (x, t ) dx dt by U x,T (dx, dt) in the last equation on p.1707 
of Meerschaert and Straka (2014). A change of variable then yields (2.7). □ 
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3. A Banach space framework 

In order to properly define the backwards and forwards equations governing the 
CTRW limits we establish a Banach space framework on which U is everywhere defined. 
Consider Co(R d x [a, b)), the space of bounded continuous functions on W l x [a, b), 
vanishing at infinity and b but not necessarily at a; i.e. the closure of the space of 
continuous functions with compact support in x [a, b) with respect to the sup norm. 
The idea is that we will consider the limit process on this space or its dual space for 
where s is the backward variable and t the forward variable. 

The crucial observation is that if f(x,cr) = 0 for all a > s, then, since D r is strictly 
increasing, 

(3.1) Uf(x, a) = 0 = T r f(x, a) 

for all a > s and r > 0. This allows us to restrict/project the semigroup {T r } r > 0 and 
all of its related operators to Co(M d x [a, b )). In particular, for / G Cq(W l x [a, b)) pick 
/ G C , 0 (lR d+1 ) such that f(x,s ) = f{x,s) for all x G R d and s G [a, b) and f(x,s) = 0 
for all s > b and all s < a — 1. Define the projection of {T r } r >o via 

T r f(x,s) := T r f (x, s) 

for all x G and s G [a, b). This is well defined by (3.1) and hence also defines a 
strongly continuous semigroup with generator A. Since U : / ha /)/° T r f dr is defined 
for any continuous function with compact support, by Fatou’s Lemma it is a bounded 
operator, and by the resolvent identity, U = —A~ x . With the same argument, T is a 
bounded operator. 

In the following we will not distinguish between T r and T r , etc. 


4. Kolmogorov backward equation 


We now define the transition kernel P for CTRW limits via 
(4-1) [ f( V )P(x, s ;dy,t) = E‘-‘[f(X,)}, 

Jy£ R d 

where / G C&(M d ). We interpret the starting point x and starting time s as the backward 
variables, and y and t as the forward variables. We also define for h(x, s ) = f(x)g(s), 

Ph(x,s):=f I P(x,s-,dy,r)h(y,T)dr= f g(T)E x ’ s [f(X r )}dT. 

Jt>s Jy£M d Jt>s 


Theorem 4.1 (Kolmogorov Backward Equation for CTRW Limits). Let h G Co(R d x 
[a, b)). Then Ph lies in the domain of A, and Ph is the unique solution to the problem 
of finding v G Cc^M^ x [a, &)) satisfying 

—Av = T/i. 


Proof. For h(x, s ) = f(x)g(s), the statement follows directly by adapting the statement 
of Proposition 2.2 onto C 0 (M d x [a, b)). For general h G C 0 (R d x [a, b)) the statement 
follows from the closedness of A, boundedness of T and the fact that functions of the 
form f(x)g(s) are a total set. Uniqueness follows from the fact that A has the bounded 
inverse —U. □ 
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Remark 4.2. Recall that in P(x, s;dy,t), we call (x,s) the “backward” variables and 
(y, t ) the “forward” variables. Unlike most backward equations, Tli 4.1 does not directly 
relate the s-derivative of the transition kernel P to the generator of spatial motion 
(acting on x). However considering the limit of solutions Ph n with h n (x, s ) = f(x)g n (s) 
and g n —> S t with supp(g n ) C (t, t + 1), by the right continuity of X t and (4.1), 

Ph n (x,s)^E x ’ s f(X t ). 

Remark 4.3. Carmi et al. (2010) derive a “backward fractional Feynman-Kac” equation, 
in the case where jumps have finite variance and are independent of the waiting times. 
In its generality, Th 4.1 above appears to be new. 


5. Fokker-Planck Equation 


In this section we show that the probability law of the CTRW limit is a unique 
solution to a FPE as long as the tail of the temporal Levy measure is time independent or 
the corresponding operator is invertible. In particular, we are interested in formulating 
the problem that is solved by the law of X t given that X s = p. 

Recall that by the Riesz Representation Theorem the dual space of Uo(R d x [a, b)) 
is the space of regular bounded measures A4(R rf x [a, b)) (Rudin, 1987) and that the 
adjoint of a densely defined linear operator A on a Banach space X is a uniquely defined 
closed operator on its dual X*. It is defined via x* G Dom(kP) if there exists y* G X* 
such that x*(Ax ) = y*(x ) for all x G Dom(A), and then A*x* = y* (Phillips, 1955). 
This is relevant as 


P*(p<8> 6 s )(dy,dt) 



P(x, s ; dy , t)y(dx) dt , 


t > s 


is the quantity of interest (its right-continuous version). 

As U and T are bounded operators, so are U* and T*. In particular, a simple 
substitution shows that 


fy*h(dy, dt) = h(dy, dt)^(y : t) + dt 


h(dy, da)H(y , a; t — a). 


' a<cr<t 


As is common, we define the convolution * in the variable t to be 


(p -k t u)(dx, dt) = 


p(x, dt — s)u(dx, ds) 


J sE[a,6) 

for every v G A4(lU x [a, b)) and family of measures {p(x, dt)} xeR d on R such that 
x i —y p(x, B) is measurable for every Borel set Bel. 


Proposition 5.1. If 7(p,t) = q(p) and H(y,t]v) = H(y,v) do not depend on t, then 
T* is one-to-one and 

^ lh = ii v *‘ h 

for h in the range of T*. The Laplace transform of the measure V ( y , •) is given by 


e V (y, dt) 


1 1 
^7 (y) + H(y,X) 
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Proof. The measures V(y,-) exist since they are renewal measures of subordinators 
with (fixed) drift 7 (y) and Levy measure h(y\ dw) (Bertoin, 1999). The statement then 
follows from basic Laplace transform theory. □ 


Theorem 5.2 (Fokker-Planck Equation for CTRW Limits). Assume T* is one-to-one. 
Let the initial condition h be given by h(dy,dt ) = p(dy)5 s (dt), where /i E A4(M d ) and 
a < s < b. Then P*h is the unique solution to the problem of finding v E A4(R d x [a, b )) 
satisfying 

= -h. 

Proof. On Co(M d x [a, b )), P<f> = f/T0 for all <fi. Hence P*h = T *U*h and equivalently, 
(ty*)~ 1 P*h = U*h for all h E A4(M d x [a, b)). Therefore (T*) _ 1 P*h is in the range of 
U* and hence (T*)” 1 P*/i E Dom(^4*) and 

A* (f$/*)~^ P* h = A*U*h = -h. 

Since A* is invertible, ^4*('L*) _1 u = 0 implies u = 0 , which implies uniqueness. □ 
Corollary 5.3. The transition kernel P(x, s;dy,t) satisfies 

—^4*(T*) _ 1 P(x, s; dy, t)dt = S x (dy)5 s (dt). 

The Fokker—Planck operator. In case that temporal and spatial jumps are uncou¬ 
pled; i.e., K is concentrated on the axes, that is 

(5.1) K{x, s, dz, dw) = K(x, s; dz x {0}) + K(x, s; {0} x dw), 

above equation simplifies further as it allows the splitting of A = T> + C into a temporal 
operator V and a spatial operator C. In particular, after integration by parts, 

d f d 

Vfix, s) = 7 (x, s)—f(x, s) + / ^-/(x, s + v)H(x, s; v) dv 

OS J v >0 OS 

and 

A fix, s) =b\x, s)d Xi f{x, s) + ^a lJ (x, s)d Xi d Xj f(x, s) 

+ I [f(x + z,s)~ f(x,s)-z l l(\\z\\<l)d Xi f(x,t)]K(x,s-,dz,{0}) 

z€M d 

Identifying Vf as \l >^f{x,s), taking adjoints we obtain 

A*f{x, t) = -^Ar*f(x, t) + £*f(x, t). 

Hence the governing equation simplifies to 

(5.2) J|-P*h = C*(^*)~ l P*h + h, 
earning C* its designation as Fokker-Planck operator. 

Remark 5.4. Under the assumption that the law of the CTRW limit has Lebesgue 
densities, (5.2) is equivalent to Equation (45) in Kolokoltsov (2009). 


The memory kernel. The non-Markovian nature of the underlying CTRW limit is 
represented by a ‘memory kernel’ as in (Sokolov and Klafter, 2006). Their Equation ( 8 ) 
corresponds to (5.2) where ('I '*) -1 “= d/dt M* t ”. This identihes the anti-derivative of 
(xj *)- 1 as the memory kernel M(t). If the coefficients of 7 (y, t) = 7 (y) and H(y, t; w) = 
H(y;w ) do not depend on t, then M = V. In many cases the measures V(y,dt) are 
Lebesgue-absolutely continuous with density v(y,t ); e.g. when 7 (y) > 0 (Bertoin, 1999, 
Prop 1.7). 


6. Anomalous Diffusion: Examples 

6.1. SubdifFusion in a time-dependent potential. Let /3 E (0,1) and define 

HAw) ' = mO)”" 3 ’ hf(w) := = fi 

We introduce the scaling parameter c > 0, and define 

(6.1) Hp(w) := 1 A [Hp(w)/c\, h c p (w) := l{w > (r(l - I3)c)~ l/I3 }hp(w)/c. 

Note that H%(w) is the tail function of a Pareto law on (0,oo), and h c ^{w) is its den¬ 
sity. This law shall be assumed for the distribution of waiting times. We also assume 
probabilities £(x,t) and r(x,t ) to jump left or right on a one-dimensional lattice. A 
CTRW with such jumps and waiting times may be represented as a Markov chain in 
M d+1 , with transition kernel 

(6.2) A c (x, s; dz, dw ) = [£(x, s + vo)8_& x {dz) + r(x, s + w)5^ x (dz )] hp(w)dw. 

Such CTRWs are a useful model for subdiffusive processes, i.e. processes whose variance 
grows slower than linearly (Metzler and Klafter, 2000). For the limit to exist asc-I 00 , 
we assume 

(6.3) £(x, s) + r(x, s) = 1, r(x, s) — £(x, s) = b(x, s)Ax. 

where b(x, s ) is a bias and Ax is the lattice spacing. The bias varies with space and 
time and is given e.g. by the concentration gradient of a chemo-attractive substance, 
which itself diffuses in space (Langlands and Henry, 2010). 

We consider the limit c —» 00 , with (Ax) 2 = 1/c. The limiting coefficients of (A r , D r ) 
are 


a(x, s ) = 1, b(x, s ) = given, y(x, s) = 0, K(x, s; dz , dw) = 5(dz)hp(w)dw, 
where 

(6.4) hp(w) = Pw~P- l l{w > 0}/r(l -13). 

and 8 denotes the Dirac measure concentrated at 0 G Apply Jacod and Shiryaev 
(2002, Th IX.4. 8 ) to see that the convergence (2.1) holds. The infinitesimal generator 
reads 

A/O, s) = b(x, s)d x f(x , s) + ^d 2 /(x, s) - dO/O, s) 
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where dl_ s f denotes the negative fractional derivative (Meerschaert and Sikorskii, 2011; 
Kolokoltsov, 2011). Given a suitable “terminal condition” / G C7(M d+1 ), the Kol¬ 
mogorov backward equation is hence 

d- s Pf (x, s ) = b(x, s)d x Pf(x, s ) + lf%Pf(x, s ) + d^_~ 1 f(x, s) 
where the negative Riemann-Liouville fractional integral of order f3 > 0 is denoted by 
(6.5) dll f{t) '■= ^ f(t + r)^- 1 dr 

(see also Bajlekova (2001)). 

For the forward equation, we note that H(x,s]w ) = Hp(w) := U! _/3 /r(l — j3) has 
Laplace transform Hp{\) = A^ -1 . Hence U(A) = A _/3 , which inverts to V(y,r ) = 
r / 3 ~ i /r(/ 5 ) — Hi_p(r). Thus (T *) -1 may be interpreted as the fractional derivative 
. The adjoint of C is given by 

£*f(dy, dt) = —d y [b(y, t)f{dy , dt)] + ^ d 2 y f(dy , dt), 


hence the distributional Fokker-Planck equation is 

d t P*[n (8) 5 S ] = —d y bd\~ h P*{(i ® 8 S ) + J^d]^ 13 P* [n <8> 5 S ) + /i (8) 


(compare Henry et al. (2010)). 


Remark 6.1. The coefficients a, b, 7 and iF above match the coefficients of the stochastic 
differential equation (7) in Magdziarz et al. (2014) where the diffusivity = 1. The 
Fokker-Planck equation also matches their equation ( 6 ). A CTRW scaling limit whose 
diffusivity varies in space and time is achieved e.g. if ( 6 . 2 ) is replaced by 

K c (x, s; dz, dw ) = A/"(gR|c -1 / 2 6(o;, s), c~ 1 a(x, s))hp(w) dw , 
where Af(dz\m, s' 2 ) denotes a univariate Gaufiian distribution with mean m and variance 


6.2. Traps of spatially varying depth. Fedotov and Falconer (2012) study CTRWs 
with spatially varying “anomalous exponent” /3(x) G (0,1). They find that in the long¬ 
time limit the (lattice) CTRW process is localized at the lattice point where /3(x) attains 
its minimum, a phenomenon termed “anomalous aggregation”. Using flux balances, 
Chechkin et al. (2005) derive a fractional diffusion equation with a “variable order” 
Riemann-Liouville derivative, which we can now rephrase in our framework. I 11 this 
example, we assume unbiased jumps of probability 1/2 to the left and right, and fix a 
Lipschitz continuous function j3{x) G (e, 1 — e) for some £ > 0. The waiting time at 
each lattice site has the density h c p, x Aw) as in (6.1), with /? replaced by /3(x). In the 
limit c —y 00 with (Ax) = 1/cwe arrive at the coefficients 

( 6 . 6 ) a(x, s) = l, b(x,s)=Q, q(x, s) = 0 , K{x\dz x dw) — 8 0 (dz)hpM(w)dw. 

As mentioned in Bass (1988, p.272), the standard Lipschitz continuity and growth 
assumptions guarantee the existence and uniqueness of a strong (pathwise) solution 
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to a stochastic differential equation with generator A given by (2.3) and (6.6). The 
negative fractional derivative of variable order j3(x) is 


d-¥ ) f(x,t) = [ [f{x,t) ~ f{x,t + w)]h/3( x )(w)dw, 

J w>0 

where h^ x ){ w ) is as in (6.4), with (3 dependent on x. As in the previous example, we 
have V(y,r) = r^^- 1 /T(f3(y)), and the Kolmogorov backward equation hence reads 

d^ ] Pf(x, s ) = i dlPf(x , s) + s) 

and the FPE 


d t P*(/a ® 8 S ) 



d] ~ m P*(jJi®S s ) 


+ n <8> S s . 


Remark 6.2. A different approach to spatially varying traps is taken in Kolokoltsov 
(2009). There, the generator 

Af(x,s) = J J i (f(x + y,s)-f(x,s))-^^S(x,s,y)d s y 

w(x,s ) f 00 . .s 1 

+ Y(^jJ 0 (f( x A + v)-f{x,s)) : ^dv 

=: Cf(x, s ) + w(x, s)d P _J(x, s ) 


for the process (A r , D r ) is assumed, where a G (0,2), j3 E (0,1 ),y = y/\y\, S(x, s, y) dsy 
is a symmetric Lebesgue-absolutely continuous measure on the unit sphere and w a mea- 
sureable function. The scaling limit process is explicitly constructed. An application 
of Theorem 4.3 therein gives 

(6.7) —d t P(y, t ) = df [w(y, t)U(y , t)] 


where U(y,t ) is the density of the potential measure of the Feller process (A r , D r ) and 
P(-,t) the probability density of X t . Assuming that w(x,s ) = w(x) does not depend 
on the time variable, we may go one step further and write the FPE for this CTRW 
limit process as 

d t P*[n(dy) ® S s (dt )] = C*P*[n(dy) ® S s (dt)\ + y(dy) <g) 6 s (dt)- 
note that, unlike in the previous example, we now have 


V(y,dt ) 


t ?-i ? 

F (P)w(y) ’ 


and C is formally self-adjoint. 


The Kolmogorov backward equation reads 

- CPf(x , s) - w(x, s)d^_ s Pf(x, s) 
where w(x, s ) may be time-dependent. 


w(x,s)d^ s 1 f(x,s), 


li 









6.3. Space- and time-dependent Levy Walks. The standard Levy Walk consists 
of i.i.d. movements with constant speed, where directions are drawn from a probability 
distribution A (d6) on the unit sphere S d_1 in M. d and movement lengths are drawn from 
a probability distribution which lies in the domain of attraction of a stable law, e.g. 
h c p(w) (6.1). We consider the case (3 E (0,1), which is termed “ballistic” since the 
second moment grows quadratically (Klafter and Sokolov, 2011). Coupled CTRWs, 
in which waiting times of length Wk come with jumps of size Jk\ = Wk, serve as an 
approximation of a Levy Walk with velocity 1. 

In this example, we consider a CTRW approximation of a Levy Walk with space- 
and time-dependent drift b(x,s). Such a CTRW is given by the Markov chain with 
transition kernel 

(6.8) K c (x : t; B x I) = / / ls(r6 l + b(x, f)/c)R(r)^(r)dr X(d9), 

Joes*- 1 Jr> o 

which converges to a limiting space-time process ( A r , D r ) with generator (2.3) and 
coefficients 

a = 0 , b l (x, s ) = given, 7 ( 2 ;, s) = 0 , 

K(x, s-BxI) = K(B x /) = [ I l s (r#)l/(r) h p {r) dr A (d9). 

Je&S 11 - 1 J r>0 

(Note that here b(x, s ) is relative to there being no cut-off function l(||z|| < 1) in (2.3).) 
The infinitesimal generator has the pseudo-differential representation (Jurlewicz et ah, 
2012; Meerschaert and Scheffler, 2008) 

Af(x,s) — b l (x,s)d Xi f(x,s)+ f f [f(x + w9,s + w) — f(x,s)]hp(w)dw\(dd) 

Joes’ 1 - 1 Jw> o 

= b\x, s)d x J(x,s) - [ (-(9, V x ) - ds) 13 f(x, s)X(d9). 

Joes’ 1 - 1 

The Kolmogorov backwards equation for the CTRW scaling limit is thus 

V(x, s)d Xi Pf(x, s) - [ (-(9, V) - d. s f Pf(x, s ) A (d9) = d^fix, s ). 

Joes' 1 - 1 

As H(x, s\ w ) = w~J/T{1 — /3) as in Example 6.1, the governing FPE is 

P*[n ® d s ](dy, dt) = —/ \i(dy ) ® S s (dt) 

The generator A does not have a decomposition into £ + V as in (5.1), and hence we 
stop here. 
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